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Abstract 

The placement of wind turbines on a given area of land such that 
the wind farm produces a maximum amount of energy is a challenging 
optimization problem. In this article, we tackle this problem, taking 
into account wake effects that are produced by the different turbines 
on the wind farm. We significantly improve upon existing results for 
the minimization of wake effects by developing a new problem-specific 
local search algorithm. One key step in the speed-up of our algorithm is 
the reduction in computation time needed to assess a given wind farm 
layout compared to previous approaches. Our new method allows the 
optimization of large real-world scenarios within a single night on a 
standard computer, whereas weeks on specialized computing servers 
were required for previous approaches. 

Keywords: wind energy, wind farm layout, optimization, opera- 
tions research 

1 Introduction 

Renewable energy plays an increasing role in the energy supply world-wide, 
and wind energy is one of the key players in the field of renewable energy. 
The wind farm layout problem entails the process of planning the placement 
of turbines (and supporting equipment) on a potential wind farm site. The 
creation of a farm layout involves the invocation of a software optimization 
module, which attempts to efficiently place the turbines while adhering to 
the constraints and optimizing the stated objectives. Often, this module 
is embedded within a specialized tool provided by wind power consultants 
such as Garrad Hassan or AWS TruePower, who offer a product such as 
Open Wind p. 

One of the problems that such tools have to deal with, apart from the 
actual optimization, is the scaling cost of computing wake effects when esti- 
mating energy capture for increasing numbers of turbines. To estimate the 
energy capture of a layout the optimization module models the free stream 
wind flowing through the site in and out of the turbines. Some degree of 



non-linear wind turbulence occurs at the outflow of a turbine and affects the 
inflow to turbines close enough behind it. Modelling this effect is necessary 
because wake has a great effect on the actual energy output of a wind farm. 
However, the time for computing the wake effect with respect to a given 
wake model such as the Park wake model [10] takes time 0(n 2 ), where n 
is the number of turbines on the wind farm. This computational effort is 
significant if one is applying iterative search algorithms such as local search, 
simulated annealing or evolutionary algorithms for the optimization of the 
placement. Such methods would need to evaluate each farm layout with 
respect to the wake model under consideration and would therefore require 
time Q(n 2 ) for each solution that is considered during the optimization pro- 
cess. 

1.1 Related Work 

The optimal siting of wind turbines on a given area of land is a complex 
optimization problem which is hard to solve by exact methods. The deci- 
sion space is non- linear with respect to how sited turbines interact, when 
considering wake loss and energy capture. Several bio-inspired computa- 
tion techniques such as evolutionary algorithms [31 0] and particle swarm 
optimization [7J have been used for the optimization. 

The different approaches for the wind farm layout problem are summa- 
rized in [11]. Wan et al, |T5l EH IH] , use cell-based approaches and compare 
different bio- inspired algorithms, each applied to the same set of wind farm 
models and parameters. They use successively more expressive layout rep- 
resentations to relax where in a cell a turbine can be located: strictly in 
the middle, anywhere, or anywhere subject to proximity constraints with 
neighbouring turbines. An alternative to cell placement was explored in [8j: 
each turbine's location is a decision variable pair of real-valued, spatial (x,y) 
coordinates. In that paper a simple evolution strategy (ES) is applied to op- 
timize the placement of the turbines. In general, an ES is effective because 
it is easily parallelized and it self-adapts the extent to which it perturbs 
decision variables when generating a new potential solution. However, the 
algorithm given in [8] is only able to deal with problems that have just a 
small number of turbines. A more powerful ES called CMA-ES has been 
used for the same problem in p3] . The advantage of this approach is that it 
allows for the effective optimization of huge layouts for up to 1000 turbines. 
Nevertheless, this approach is still computationally expensive, requiring up 
to weeks to solve for large wind farm layouts on computing servers. 

Many of the above-mentioned approaches make simplified assumptions 
regarding the realism of the wake models and the wind scenarios. Even 
in recent work, it is either completely ignored, assumed to be constant, or 
highly simplified wind scenarios are used (except for [HUH]). For example, 
in the optimization formulation in [9], the turbines are placed in regular 
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grids, and wake is considered by enforcing minimal and maximal distances 
between the turbines, thus effectively neglecting physical effects. In [13] . 
distance-dependent wakes are considered, but the underlying wind scenario 
was randomly generated, and is the result of unrealistic assumptions. On 
the other hand, the above-mentioned industry tool AWS OpenWind contains 
elaborate wake models that are used for real-world scenarios, but its simple 
optimizer neither considers the turbines' vicinities when perturbing layouts, 
nor adapts parameters for the perturbation during the run. 

Note that, placing turbines on a defined area is loosely related to the 
difficult problem of packing discs in shapes [U [12]. There, the task is to ar- 
range n identical discs without overlap entirely inside a square on the plane, 
such that the discs have the largest possible diameter. As the turbines influ- 
ence each other via non-linear wake-effects, these "influence areas" can be 
assumed to be circular for very simple scenarios of uniform wind distribu- 
tions^] The naive goal then can be to minimize the overlap of the influence 
discs. However, the theoretical results from the disc packing are not appli- 
cable in our context, as it would be extremely simplifying to assume that 
these discs are circular due to non- uniform wind distributions. Addition- 
ally, the sizes and irregular shapes are highly dependent on the interaction 
of the turbines: due to interactions, the wind resources available at a turbine 
change, and so does then the influence area around a turbine. 

1.2 Our Contribution 

In this paper we describe the design of a fast and effective randomized local 
search algorithm to act as a farm layout optimization module. The algorithm 

1. uses a fast approach for evaluating new layouts, 

2. makes effective use of the problem characteristics to produce new lay- 
outs, and 

3. can efficiently handle infeasible regions in order to respect geographical 
constraints. 

To be precise, where previous algorithms required time fi(n 2 ) for the 
computation of the wake effects within a new layout, our algorithm requires 
@(kn), where n is the number of turbines on the wind farm, and k the 
number of turbines whose location was changed over the previous layout. 

The new local search algorithm is capable of generating good layouts for 
hundreds of turbines quickly on standard computers, without the need of 
specialised computing hardware. This allows it to accommodate the emerg- 
ing requirements of larger farms. For example, the Horse Hollow Wind 

x Note that these discs are not necessarily of infinite size, as a cut-off distance can be 
set once the influence becomes negligible. 
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Algorithm 1: Procedure for evaluation of wake effects due to the Park 
model |8j 

1 Given {X, Y} as turbine locations, turbine thrust coefficient Ct, rotor 
diameter R, landscape-specific wake spreading factor k ; 

2 a = 1 — v 7 ! — Ct, b = k/ R, u unit step function, 
o = {xi — xj)cos9 + (yi — yj)sin9; 

3 d 



\o\\, a = tan k 



4 for i = 1 to number of turbines do 

5 for 6 = 0° to 360° do 

6 for j = 1 to n-1 and j ^ i do 



9 
10 



cos 



o+R/k 



I (xi—Xj + ^cos9) 2 +(yi—yj + ^sin6 



_ Vdef {i}j) = u(5ij - a) {1+b a di ,y 2 ; 
Vdeff - 



2. 



d(0) = *(6) x (1 - Vdefi); 



Energy Center in Texas, USA operates with 735.5 megawatt capacity and 
consists of more than 300 turbines spread over nearly 47,000 acres (190 
km 2 ). We compare our new algorithm to two state-of-the-art approaches: 
the CMA-ES approach presented in p3] and the local search approach used 
in the industry tool AWS OpenWind. Our experiments show that the new 
algorithm outperforms both approaches in terms of the quality of the results 
and the running time. 

We proceed as follows: Section [2] provides a description of the farm 
layout optimization problem. For this, it is necessary to explain how wake 
affects wind characteristics as wind propagates through the site, and how 
the expected energy capture is calculated. Section [3] describes our algorithm. 
Section [4] gives details about our experiments, and discusses the outcomes. 
Then, the capability to deal with infeasible areas is shown, using a scenario 
based on an existing near-shore wind farm. Finally, Section [6] summarizes 
our findings. 

2 Wind Farm Layout Optimization 

Let X = {xi, . . . , x n } and Y = {y±, . . . , y n } be a set of x and y coordinates 
of n wind turbines in the plane. The goal is to find a set of coordinates such 
that the energy output of the whole wind farm is maximized. Furthermore, 
the layout has to respect certain constraints. 
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2.1 Energy Output 

Depending on the chosen coordinates the overall energy output of the wind 
farm varies as we have to take the wake effects into account. 

We consider the Park wake model. In this model the wake effects on a 
turbine i change the wind resource available to it along different directions 
by reducing the scale parameter c of the Weibull distribution estimated for 
the entire farm, which is also called the freestream wind resource. This is 
dependent on its location and the location of the rest of the turbines. Hence, 
we have a parameter q for each turbine i: its computation is complex and 
involves wind velocity deficits Vdefj that the turbine i experiences due to 
the influence of other turbines j, j ^ i (see Algorithm [I]) . We refer the 
reader to [8J for a detailed presentation on the computation of this parameter 
when considering wake effects in the Park wake model. The expected energy 
output 7) of the whole wind farm is given by 

E farm [7]] = £ ^ P(0) jf p(t;(0), ^ Y), k(9))(3\v). (1) 

In this equation v is the wind speed, and the function f3 l (v) defines the 
power curve for turbine i. Wind speed v however is a random variable with a 
Weibull distribution, p(v(9),Ci(9, X, Y),k(9)), which is estimated from wind 
resource data and considers the wake effect using X and Y. This distribution 
is also a function of the wind direction, 9 which varies from 0° -360°. Note 
that the shape parameter of the Weibull distribution is not influenced by the 
Wake effects here. Additionally, wind flows from a certain direction with 
some probability P{9). 

2.2 Constraints and Assumptions 

We have the following constraints placed on our optimization function. The 
first one enforces an upper bound on the area of the farm. This constraint 
ensures that we can only place a turbine i within a certain area, which is a 
realistic constraint for most layout problems. For a rectangular farm with 
length I and width w this constraint is satisfied iff < x\ < I & < yi < 
w,Vi. 

The second constraint regulates the spatial proximity, as it dictates the 
minimal distance within which two turbines can be set up. It is satisfied iff 
yj (xi — Xj) 2 + (yi — yj) 2 > MR, Mi Vj where R is the rotor radius and M is 
a proximity factor usually decided ahead of the optimization based on the 
make and model of the turbines used. We use M = 8 based on the industry 
standard. 

In addition to the above constraints, we assume that all turbines have the 
same power curves (approximated as piecewise linear functions) and that the 
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same wind resource spans the entire farmj^] The assumptions can be revised 
in a very straight forward manner to generate more realistic scenarios. 



3 Turbine Distribution Algorithm 

In this section, we describe our algorithm for the optimization of wind farm 
layouts. It includes a problem specific local search operator and a faster 
evaluation. 



3.1 Computational Speedup 

First, we describe how computation time can be saved, if the new layout 
differs from the old one only in the location of a single turbine. 

For the computation of the (turbine-specific) wind resources, the Weibull 
scale parameters are adjusted up to 24 x n 2 times (line 10 of Algorithm [l]), as 
the wind direction is discretized into 24 wind directions in our model, and the 
mutual influence of the n turbines has to be considered. The computation 
of the influence matrix results in an evaluation time that is quadratic with 
the number of turbines, causing evaluation times of up to 30 seconds on 
standard hardware for a 1000 turbine layoutj^] 

As our optimization algorithm modifies only a single turbine of the cur- 
rent layout, the evaluation can be speeded up by updating the velocity 
deficits in line 10 of Algorithm [T] intelligently. First, for a moved turbine i 
the influence of all other turbines on i has to be computed conventionally, 
requiring n-1 influence checks. Second, the other turbines' velocity deficits 
have to be updated as i's influence on these may have changed. This can be 
done as follows. For each such unmoved turbine j: 



T r, t 9,NEW /7t/j f 8,OLD\ 2 / T r, f e,OLD\ 2 . ( ~\r j r9,NEW\ 2 fo x 

Vdefy =^[Vdef- ) -{Vdeffa ) + [Vdeffa ) (2) 

where Vdef e - ,OLD is the velocity deficit the turbine j experiences in the old 

layout, and Vdef®j^ LD is the influence that i had on j, for a given wind 
direction 9. 

Over the course of the algorithm's run, we only have to initially create the 
three dimensional array carrying all mutual influences for all wind directions 
once. In subsequent evaluations, we can use and update this matrix to 
speedup the evaluations. Generalized, this allows for a resulting runtime for 
such a subsequent evaluation of only Q(kn), when the locations of k out of 
n turbines were changed]^] 



To increase accuracy, these resources can be estimated for different parts in the farm. 
3 This test and all subsequent experiments were performed on AMD Opteron 250 CPUs 
(2.4GHz), on Debian GNU/Linux 5.0.8, with Java SE RE 1.6. 

4 We use k — 1 in our experiments, which results in a significant speedup (see TableJTJ). 
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3.2 The Algorithm 

It is clear that the constraints discussed in Section 12.21 are vital to the con- 
struction of the algorithm. To ensure constraint handling in the optimiza- 
tion, a random local search algorithm was purpose-built for the application. 
The turbine distribution algorithm (TDA) described in Algorithm [2] itera- 
tively displaces a single turbine in order to increase the energy gain while 
ensuring constraints are upheld. 

In order to ensure that the turbines' initial placement respects the safety 
distance constraint, our TDA follows the approach taken in |14j . There, the 
algorithm deterministically initializes the turbines in a grid formation of 
greatest space. The grid is constructed in such a way that the distance 
between the columns and rows is maximized, including the placement of 
turbines on the borders of the wind farm area. This is a straightforward 
approach, which is used frequently in practice, as the wake effect is already 
reduced to some extent (when compared to tighter layouts), even without 
considering the directional distribution of the wind. 

Over the course of the optimization, new layouts are created based on 
the best-so-far turbine configuration. For the next layout, a copy of the best- 
so-far is made, and a modification is then applied, in which a single turbine 
is selected uniformly at random and is shifted by a displacement vector if. 
In the following, we motivate our way to compute this displacement vector. 

Initially, the vector v is determined by the normalised sum of the differ- 
ence in distance between the current turbine and its nn nearest neighbours. 
The displacement of the current turbine by v would move the turbine away 
from the nearest neighbours. In order to ensure non-deterministic perfor- 
mance within sensible limits, TDA applies a normal distribution with a 
fixed standard deviation to the direction of the vector and uses an adapt- 
able standard deviation on a mean over a normal distribution to determine 
the distance for the vector. As a superior energy output for the wind farm 
could be obtained by moving turbines closer together, this new vector may 
be randomly reversed with a preset probability. 

In the next step, the new location of the turbine after displacement by 
the resultant vector v 1 is investigated to ensure that it would not be placed in 
an illegal position, currently defined as outside the wind farm bounds or too 
close to another turbine. If the turbine would be in an illegal position, the 
distance of v 1 is reduced to a legal value. An example of this displacement 
using nn = 2 may be seen in Figure [TJ Note that, in this example, the 
two nearest neighbours of the current turbine are turbines 3 and 4, and 
hence these are the only turbine locations that affect the new position of 
the current turbine. 

At the conclusion of this layout modification, the quality of the new 
layout is compared against that of the best-so- far, using Equation [j] and the 
computational speedup presented above. If the new potential solution is of 
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Algorithm 2: Turbine Distribution Algorithm (TDA) 

1 Given number of turbines n, map bounds x max , y m ax, nearest 
neighbours nn, reversal probability p, displacement distance standard 
deviation (J[dis} k for k G {1, . . . , n}, and direction standard deviation 
a {dir} k for k e {1, ... , n}, place 

{X, Y} = {[0, . . . , x ma J, [0, . . . , y max ]} n in the grid formation of 
greatest space; 

2 Determine f({X,Y}) as per Equation [IJ 

3 for i = 1 to number of evaluations do 
Set {X,Y}' = {X,Y}; 

Select turbine k uniformly at random for the i th modification, 
denote as {x' k , y' k }; 

Determine nearest neighbours: {x'{, y'{}, {x" n , y" n }; 

2^j=i x k ~ 00 j ' 2^j=i Vk ~ Vj 
v = v/\\v\\; 

Select 9 normally distributed over fi = Zv and a = cr^ r y , ; 
Select d normally distributed over \i = and a = cr^ dis y k ; 
if = 6*d; 

Set v , =—v i with probability p, as tighter groups may increase the 
farm's output; 

If applying if would place {x' k ,y' k } in illegal area, reduce length of 
if until legal; 
Displace {x' k ,y' k } by if; 

Determine f({X,Y}') as per Equations [l] and [2J 
if f({X,Y})<f({X,Y}>) then 
Set {X,Y} = {X,Y}'; 
Increase cr{ dis } fc ; 

else 

_ Decrease cr {dis}k ; 
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Figure 1: Illustration of vector displacement using TDA with nn = 2. 

higher quality, it replaces the previous best. 

Furthermore, when the quality of the layout improves due to the dis- 
placement of a turbine, the above-mentioned adaptable standard deviation 
for the distance of that particular turbine increases to allow for increased 
exploration. Similarly, when the new layout is of lower quality due to the 
displacement of a turbine, the standard deviation of this particular turbine 
decreases to allow displacement exploitation. Effectively, this gives TDA 
the ability to autonomously switch between exploration and exploitation. 
There will be local optima due to the local search properties of the algo- 
rithm, and our strategy to increase the potential of escape is the combination 
of n adaptable displacement parameters, in combination with values drawn 
randomly from different normal distributions. 

4 Experimental Investigations 

In order to justify our design decisions, we performed experimental investi- 
gations on which we report in the following. First, we introduce the scenario 
that defines the wind resource present at an imaginary prospective site of 
a wind farm. Then, we describe TDA's parameter settings and the test 
scenarios for the subsequent comparative study. There, we compare the so- 
lution quality and runtime of our algorithm to a tuned version from |14J and 
the industry tool AWS OpenWind [I] . 

4.1 Algorithm Settings and Scenario Information 

To evaluate our algorithm's performance we set up a realistic scenario with 
the wind resources defined as Scenario 2 in [H] . The prevailing wind direction 
covers a broad sector of about 105°, and the wind intensity per direction 
is given by Weibull distributions. This results in non-zero probabilities for 
wind coming from any direction, and therefore, one has to optimize the 
layout to work with minimum wake loss along all the wind directions. 

For the subsequent experiments, the internal parameters of our algorithm 
were set as follows. 
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Figure 2: Study: influence of the number of considered nearest neighbours 
on the performance. Shown is the energy gain and its standard deviation 
achieved over the initial layout. 



The direction standard deviation is set to cr^ir}- = tt/6 for all i tur- 
bines. Thus, the resulting displacement direction will be within ±30° with a 
probability of 68.2% of all directional adjustments, and within ±60° with a 
probability of 95.4%. Furthermore, the reversal probability of the displace- 
ment was set to p = 0.2. 

The displacement distance standard deviation (J{dis} i ls se t and adjusted 
by the algorithm for each turbine individually, depending on the initial lay- 
out. If the minimal initial distance between turbines is d, then cr^ dis y. = 
(d — 8i?)/3. Thus, we restrict the displacements to small distances, while 
considering the safety distance of 8R between the turbines. 

As described in Section 3J2 these parameters are adjusted over the course 
of the optimization depending on whether recent displacements were success- 
ful or not. This allows for an efficient exploration as well as exploitation. 



4.2 Impact of the Nearest Turbines 

In order to understand the influence of the number of nearest neighbouring 
turbines nn that influence the relocation of a chosen turbine, we ran an 
initial set of experiments. We chose the scenarios with 40, 100, and 400 
turbines, and ran each experiment with a budget of 10,000 evaluations, nn 
varied from 1 to 8, and the results of 100 repetitions per scenario are shown 
in Figure [2} 

We noticed that the number of nearest neighbours sometimes has some 
influence on the performance of the algorithm. Depending on the scenario, 
a smaller value of nn seems to be beneficial. For example for the 100 tur- 
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bines scenarios, using nn = 2 instead of nn = 8 yielded in layouts with an 
increased energy output of 1.1% on average. For the 40 turbine scenario, 
however, the maximum distance in the averages is only 0.6%, and for the 
400 turbines just 0.1%. A possible explanation is that with four and eight 
neighbours, the influence of two opposing neighbouring turbines is cancelled 
out: the resulting displacement vector for a selected turbine just has a 'un- 
specific' direction, instead of an 'explorative' direction as is the case when 
only one or two neighbours are considered. 



4.3 Experimental Results 

To assess the performance of our algorithm, we compare it directly with the 
results from [J4] and the approach implemented in AWS Open Wind. In our 
computational study, we set up several scenarios with varying numbers of 
turbines, and varying farm sizes. For n = 10, 20, . . . , 100 a quadratic farm of 
3 x 3km 2 was chosen, and for n = 200, 300, 400, 500, 1000 rectangular farms 
of 8 x 5km 2 , 10 x 6km 2 , 12 x 6km 2 , 14 x 7km 2 , and 20 x 10km 2 . 

To allow for a better comparison with the CMA-ES algorithm presented 
in |14j . where 10,000 generations with a total of 200,000 evaluations were 
performed, we ran our algorithm for 10,000 and 200,000 evaluations — we 
stopped the algorithms if no improvement was reached after 1,000 evalu- 
ations. In order to compare our TDA not only to an academic system, 
we additionally implemented the optimizer component that is part of the 
industrial tool AWS OpenWind, and ran it for a maximum of 200,000 eval- 
uations. Quickly, we noticed that given the wide range of scenarios, the 
original CMA-ES fails to produce consistent results — in particular, it per- 
formed poorly for 70 and 80 turbines. Despite the fact that it has elaborate 
adjustment capabilities, it is not always capable of learning the problem- 
specific features of the solution space within reasonable run-times. Through 
the increase of CMA-ES's initial standard deviation, the explorative phase 
was extended, and we were able to increase its quality performance on these 
two particular scenarios to a level comparable with TDA. In the following, 
we refer to this tuned version as CMA-ES*. 

The initial layout, where turbines are placed in a maximally spaced grid, 
was identical for all algorithms. Each scenario was repeated 30 times, and 
the results are listed in Table [l] and in Figure [3p1 As can be seen, our 



algorithm's results (when using only 10,000 evaluations) are comparable 
to those of CMA-ES*, and in many cases have a higher solution quality. 
Additionally, our algorithm just uses a fraction of the evaluations and of the 
time to reach the results. Exemplarily, the academic CMA-ES* needs 2h 
for the 40 turbines scenario, while even a better gain is reached within just 
three minutes by our algorithm. And for the 400 turbines, our algorithm 



5 In the scenarios with 10 and 20 turbines, the turbines are very far away from each 
other, with just minimal wake effects occurring. 



11 



Table 1: Results. Reported are the predicted energy outputs in kW (average 
and maximum values, and the standard deviation), time in hours. The 
algorithm with the best average energy output has been highlighted in bold, 
the second-best in italics. 



n Initial 



CMA-ES* 200k 
avg s tdcv max time 



Open Wind 200k 



avg s tdcv 



max time 



10 7.290E+4 

20 1.448E+5 

30 2.015E+5 

40 2.630E+5 

50 3.247E+5 

60 3.688E+5 

70 4.341E+5 

80 4.707E+5 

90 5.207E+5 

100 5.535E+5 

200 1.248E+6 

300 1.861E+6 

400 2.415E+6 

500 3.062E+6 

1000 6.023E+6 



7.306E+4 3 . 63E+1 
1.448E+5 .oos+o 

2.096E+5i.73B+2 
2.742E+5 2 .07B+2 
3.367E+5 3 .01B+2 
3.947E+5 3 .47i3+2 
4.435E+5 8 . 98i3+ 3 
5.077E+5 imE+3 
5.541E+5i .22B+3 
5.994E+5 2 .i 8£;+ 3 

1.301E+6 4 .36B+2 
1.935E+6 7 . 76£ ; +2 

2.549E+6i. 2B+3 
3.196E+6 9 .63+B3 
6.298E+6 2 . 36B+4 



7.315E+4 0.10 
1.448E+5 0.05 
2.100E+5 1.29 
2.746E+5 2.01 
3.372E+5 3.82 
3.953E+5 5.45 
4.559E+5 5.90 
5.098E+5 9.78 
5.567E+5 9.64 
6.029E+5 14.0 
1.302E+6 40.1 
1.937E+6 111 
2.550E+6 221 
3.201E+6 324 
6.335E+6 327 



7.290E+4 .oo£+o 
1.448E+5o.oob+o 

2.052E+5i.i7B+3 
2.671E+5 8 .08E+2 
3.269E+5 7 .28B+2 
3.837E+5i.i9B+3 
4.370E+5 9 . 03S+ 2 
4.884E+5i.i8£,+3 
5.333E+5i.o2E+3 
5.762E+5i.48B+3 
1.273E+6 2 .63E+3 

1.893E+62.77B+3 

2.480E+6 2 . 2 oB+3 
3.H8E+63.55B+3 
6.202E+6 4 .8i B+ 3 



7.290E+4 0.01 
1.448E+5 0.01 
2.067E+5 0.01 
2.688E+5 0.02 
3.282E+5 0.03 
3.858E+5 0.04 
4.389E+5 0.07 
4.906E+5 0.08 
5.356E+5 0.10 
5.784E+5 0.14 
1.278E+6 0.39 
1.898E+6 0.87 
2.485E+6 1.46 
3.124E+6 2.26 
6.210E+6 9.43 



TDA 10k 

avgstdcv max time 



avg s tdev 



TDA 200k 
max time 



loss 



gain 



10 

20 

30 

40 

50 

60 

70 

80 

90 
100 
200 
300 
400 
500 
1000 



7.305E+4L73B+1 
1.448E+5 7 .o5E-i 

2.123E+5 3 . 74E +2 

2.772E+5 4AeE +2 
3.392E+5 4mE+2 
3.980E+5 b .25E+2 
4.512E+5i. WE +3 
5.044E+5i.i6s+3 

5.548E+5 1.34E+3 
6.015E+5LS2E+3 

1.309E+6 S .99S+2 
1.949E+6 1 . 37E+3 

2.55SE+6L62E+3 

3.211E+6 imE+3 

6.363E+6 2 A9E+3 



7.308E+4 0.01 7. 
1.448E+5 0.02 1. 
2.131E+5 0.04 2. 
2.781E+5 0.05 2. 
3.401E+5 0.06 3. 
3.991E+5 0.08 4. 
4.537E+5 0.09 4. 
5.083E+5 0.11 5. 
5.569E+5 0.12 5. 
6.041E+5 0.14 6. 
1.311E+6 0.33 1. 
1.952E+6 0.55 1. 
2.557E+6 0.84 2. 
3.215E+6 1.19 3. 
6.368E+6 3.69 6. 



309E+4 2 . 47£;+1 7. 
448E+5i.3ob+i 1. 
135E+54.08B+2 2. 
791E+5 5 .75B+2 2. 

412E+5 3 . 42 £;+2 3. 
011E+5 5 .25i5+2 4. 

555E+5i.57£+3 4. 
090E+6i.oib+3 5. 
609E+5 9 .38i5+2 5. 
083E+5i.28£;+3 6. 
323E+6 9 .16B+2 J- 

971E+6 9 . 72 £;+2 1 

584E+6i.i2£+3 2. 
249E+6i.55J5+3 3. 
,449E+6 2 .04i3+3 6. 



314E+4 0.22 
449E+5 0.50 
144E+5 0.75 
806E+5 1.02 
418E+5 1.29 
022E+5 1.58 
591E+5 1.86 
108E+5 2.15 
625E+5 2.52 
113E+5 2.83 
325E+6 6.73 
973E+6 11.3 
586E+6 17.0 
251E+6 24.5 
454E+6 75.0 



0.1% 

I. 0% 
2.7% 
4.6% 
6.7% 
8.6% 

II. 1% 
13.0% 
14.8% 
16.9% 
9.6% 
10.2% 
11.7% 
11.2% 
11.9% 



100.0% 
98.1% 
92.2% 
88.7% 
84.8% 
80.4% 
72.8% 
72.8% 
68.9% 
63.9% 
96.8% 
87.4% 
82.6% 
90.1% 
86.2% 
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Figure 3: Performance comparison: our algorithm TDA (using 10,000 and 
200,000 evaluations) versus the tuned academic CMA-ES* and the industry 
tool OpenWind (both using 200,000 evaluations). Shown is the energy gain 
achieved over the initial layout, and the standard deviation over the 30 
repetitions. 

needs less than an hour to achieve the results for which CMA-ES* requires 
> 200/i. If given 20 times the number of evaluations, our algorithm is able 
to produce layouts that produce up to 1.4% more energy than the CMA-ES* 
based algorithm from [14J. In addition, it is significantly faster. 

Independent of the scenario, Open Wind's optimizer achieves an average 
improvement over the initial layouts of only 2%. Responsible for this is 
most likely the lack of self-adaptation, as the optimizer gets stuck in local 
optima. This happens very quickly: typically within the first 1,000 layout 
assessments, which is reflected in the short running times. Contrary to this, 
our TDA adapts the relocation parameters for each turbine independently 
of the others. The effects of this ability to change from exploration to 
exploitation and back is reflected in the large performance advantage of 
TDA over OpenWind. 

To the results, we added some derived statistics. Pi oss = — wake los , s , 

' luis energy captured 

denotes the average percentage of unused wind energy due to wake effects. 
As expected, this loss increases for the scenarios with 10-100 turbines, as 
more turbines are packed into the same area, inflicting increased wake losses. 
For 200-1000 turbines, this value is relatively constant, as the chosen sizes 
of the farms vary such that the farms could contain roughly 2n turbines. 
Pgain is the average percentage in energy gained over the previous scenario 
(i.e. after increasing the number of turbines by 10, 100, or 500). Again, the 
effect of the increased mutual wakes is reflected in a decreasing benefit of 
adding turbines, when compared to the next smaller scenario. 
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Figure 4: Demonstration of the TDA infeasible area modelling capability. 
The image on the left portrays a satellite image of the Woolnorth wind 
farm in Tasmania, Australia [S]. The right images are examples of the the 
loose adaptation used in the modelling tool, and from left to right model 
the scenario at (252.9 MW), 5,000 (264.91 MW), and 20,000 (265.8 MW) 
evaluations. 

5 A Real- World Problem: Dealing With Infeasi- 
ble Areas 

It must be noted that the layout of wind farms are generally unable to 
be denoted simply by a rectangular area specified only by a width and 
height. Actual wind farms may have uneven boundaries due to proximity to 
unstable ground, lakes, cliffs, or may simply not have authority or ownership 
to build in certain locations. These and other geographical constraints may 
additionally exist within the bounds of the wind farm. Any algorithm that 
attempts to realistically model a wind farm must model these constraints. 
TDA has been adapted to model these infeasible areas, adding an extra 



constraint to those specified in Section 2.2 By modelling these infeasible 
regions as sets of shapes, any operation which would move a turbine into 
an infeasible region is corrected in the same manner in which the proximity 
and farm area constraints are handled above: when the application of the 
displacement vector v 1 (line 13 of Algorithm [2]) would place the turbine in 
an infeasible zone, the distance of v 1 is reduced until the displacement is 
legal. 

This extension has resulted in a capable wind farm modelling tool. In 
order to demonstrate the ability of the algorithm to model real-world wind 
farms, Figure [4] presents both the satellite image and a loosely modelled 
representation of the Woolnorth wind farm in Tasmania, Australia (40.685°S 
144.717°E). The modelled scenario uses the 37 turbines contained in the 
north Woolnorth site. However this scenario is merely a proof of concept, 
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as it is not using the wind characteristics of Woolnorth, nor specific internal 
map or terrain information of the actual site. Only the coastal details have 
been represented, and the above-mentioned wind resource is used (Scenario 
2 in 0). 

As can be seen, the turbines observe the constraints of the problem 
and are quickly distributed into a superior formation from evaluation 
through to evaluation 20,000. As the wind is predominantly from the west- 
ern direction (between 120° to 225°), the turbines tend to form in staggered 
north/south columns while leaving space along the east/west directions. 



6 Conclusion 

In this paper, we have presented a fast and efficient algorithm for the lay- 
out optimization of large wind farms. It takes problem-specific features into 
account, and benefits from the achieved reduced computational complexity 
of a layout evaluation when considering the Park wake model. As a re- 
sult, our algorithm achieves higher quality results than existing approaches, 
while the assessment speed-up allows for an optimization within minutes or 
hours instead of days or weeks (effective speed-up factors of up to 270 were 
observed). As the parallelization of the neighbourhood investigations is nat- 
ural with most local search approaches, we expect that a parallelization of 
our approach can yield further speedups. 

Although we considered one specific wake model, namely the Park wake 
model, it is important to note that our optimization algorithm can be easily 
applied to other wake models such as the deep array wake model [2]. 
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